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Breached pairing superfluidity: Possible realization in QCD 
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We propose a wide universality class of gapless superfluids, and analyze a limit that might be 
realized in quark matter at intermediate densities. In the breached pairing color superconducting 
phase heavy s-quarks, with a small Fermi surface, pair with light u or d quarks. The groundstate 
has a superfiuid and a normal Fermi component simultaneously. We expect a second order phase 
transition, as a function of increasing density, from the breached pairing phase to the conventional 
color-flavor locked (CFL) phase. 
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Because the primary one-gluon exchange interaction 
between high-momentum quarks is attractive for quarks 
in the color antisymmetric 3 channel, it is a firm pre- 
diction of QCD that cold dense quark matter is a color 
superconductor. At asymptotic densities (ignoring the 
c, b and t quarks) the ground state is well understood: 
quarks of all three flavors it, d, and s, pair according 
to the BCS mechanism, forming the color-flavor locked 
(CFL) phase 0]. 

It is much less clear what QCD predicts for the ground 
state at subasymptotic densities, which could be relevant 
for describing neutron stars. Differences among the quark 
masses cause mismatches among the Fermi surfaces of 
the species which potentially pair. There is no longer an 
abundance of degenerate low-energy particle-particle or 
hole-hole states with opposite momenta both near their 
Fermi surfaces, and so it is less obvious what modes are 
the best candidates for coherent alignment by attractive 
interactions. 

One much-discussed possibility is the LOFF phase 
(Larkin-Ovchinnikov-Fulde-Ferrel 2j). In the context of 
QCD these ideas lead to the crystalline color supercon- 
ductivity Here we suggest quite a different possibility 
for ordering with mismatched Fermi surfaces, that might 
be realized at intermediate densities in QCD. 

In QCD at intermediate densities, /i ~ 200 — 300 MeV, 
there is not only mismatch in quark Fermi surfaces but 
also a different dispersion relation, since the heavy s- 
quark, unlike the light u and d quarks, needs not be 
ultra-relativistic. This difference makes plausible a pair- 
ing phase, wherein strange quarks are raised to higher 
kinetic energies to exploit favorable possibilities for cor- 
relation energy through pairing. The opposite limit - 
pairing of a heavy and light species when heavy species 
has the larger Fermi surface — is the arena for the inte- 
rior gap phase discussed in 

For illustrative purposes we analyze a toy model with 
a massive s and a massless it quark. The Fermi momenta 
are related to chemical potentials as 



p F =p,- 6p, e 



Pf 



V(m + M 2 -^ (i) 



where Sp e will be tuned to enforce number equality (a 
stand-in for electric neutrality) . We are interested in the 



case when the Fermi momentum for the s-quark is smaller 
than that for the it-quark, p s F < p^. For simplicity, we 
shall linearize the s quark dispersion near its Fermi sur- 
face. We have checked that this simplification does not 
alter our results qualitatively. Our simplified model then 
has dispersion relations 



V u (p-p F ), 



v s (p- Pf ) 



(2) 



with V s < V u . 

In promoting particles of the heavy species to pair 
around the large Fermi surface of light species, there 
are two competing energetic factors to consider. These 
are the single-paticle energy cost of such promotion, 
V s \p F — p F \ per pair, versus the gain from creating a 
pair, k(V w + V s ), where k is the momentum gap. There 
is a net profit when 
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For flat dispersion for the s-quark, V s -C V u , promotion 
of s-quarks to a higher it-quark Fermi surface does not 
cost much energy, and a paired phase is favored. Trans- 
lating this into possible effects of non-zero s-quark mass 
the allowed range for pairing is given by 
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It can be more economical to promote heavy parti- 
cles to higher Fermi momentum than to equalize the two 
Fermi surfaces by deforming both. In earlier work Q 
we suggested the terminology "interior gap" for this phe- 
nomenon, motivated by the peaking of the gap parame- 
ter interior to the large Fermi surface that carries most 
of the spectral weight for gapless modes. Following that 
thought, it would be natural to call the limit we are an- 
alyzing here "exterior gap" . But to emphasize the com- 
mon physical mechanism, we use "breached pairing" to 
cover both. The pairing is breached, in that it vanishes 
in a solid annulus in momentum space. 

To allow quantitative estimates in analytic form we 
consider a schematic model involving pairing of u and s 
quarks. Starting with non- interacting degenerate Fermi 
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gases of u and s-quarks, turn on a weak attractive inter- 
action between light and heavy species with a coupling 
— g < 0. In a basis of light particles and heavy holes the 
quadratic part of the Hamiltonian is 



Hquad = ( 1pt p Ipa-p) 



"A -E* 



s-p 



(5) 



where the gap parameter is defined as A — 
g(2ir)~ 3 J d 3 p(t[)'l p ipl_ p ) bp, where the groundstate is as- 
sumed to be the breached pairing state. W qua d can be 
diagonalized by the Bogoliubov transformation 



tpup 
rs-p 



sm( 



cos f 



V y s — p 



(6) 



A/ 



A 2 . The new fermion fields ip 



with sin 20 F 

define the form of the interior gap wavefunction up to a 
relative phase (due to degeneracy) by i/j\0)bp = 0. In the 
meantime, we obtain two branches of quasi-particle exci- 

e D ± Jet 2 + A 2 , where 



tations with the spectra E — o p _i_ y Cp 

± _ i 



e 



p p 

^(e p ±e p ). The order parameter is related to 9 p by 



{^up^l-p) bp = \ sin2# p outside the breached region. 

The energy of the breached pairing state \0)bp can be 
computed from the diagonalized Hamiltonian as 



(H) BP = -A 1 
9 



D 



d 3 p 

(2tt) 3 



.+ 2 



(7) 



The first term is the mean field potential and the inte- 
gration is restricted to the area D defined by (p — A) < 
|p| < Pa an d Pa — IpI — (Po + where p^ are the two 
roots of — 0, po is the momentum defined by e 



to = * 



Pa=P F ~ Spf 



V s 



yu + ys 



(8) 



and A is the ultraviolet cutoff. (We will always be assum- 
ing po > A > A, Sp F .) Thus 



P A =p F -^p F ±-dSp F - — 



(9) 



with SpF = pp—pp > 0. Varying the groundstate energy 
with respect to A, d(H) bp /dA — 0, we find the integral 
equation for the gap parameter 



1 = .9 



d 3 p 1 



(10) 



Note that the integrand peaks near to |p| = po where the 
energy denominator vanishes for A — > 0. 

We can illuminate the physical meaning of p^ . Pairing 
of u and s quarks produces two quasiparticle excitation 
branches. While E~ < always holds, the E+ branch 
changes sign in the region [pXjPaI- The negativity of 



E+ means that the corresponding states are singly occu- 
pied by, in our chosen particle-hole basis, the u-quarks. 
Since this region has single occupation, it does not con- 
tribute to the pairing interaction, and it must therefore 
be excluded from the gap equation integral. The same 
phenomenon appears in a more formal treatment with 
Green functions |5j . Integration over energy yields a non- 
zero contribution only if the two poles of the Bogoliubov 
quasiparticles reside in different half planes of the whole 
complex energy plane. 

Two cases can arise, depending on whether the E+ 
branch can intersects the zero energy axis: 



BCS superfluidity. For 5p F < 



7#F' E p never in " 
tersects E+ = 0. Then E+ > and E p < for all p. 

Pairing is possible in the whole range (po — A) < |p| < 

(po + A). We arrive at the BCS solution for the gap 



A = X(V U + V s ) exp - 



N+(0)g 



(11) 



where A^ + (0) = / d 3 p<5(e+) is the density of states, 

and A is the ultraviolet cutoff. Since N + (0) ~ p 2 ,, the 
gap parameter decreases as a function of increasing Spp 
with p F fixed. 

Breached pairing (BP) superfluidity. This occurs 
when SpF > ^y^ Vs ■ E p = has two roots. We in- 
tegrate below p^ and above p^ , but exclude the region 
Pa — IpI < Pa- A solution of the gap equation exists 
only if g > g c , with the critical coupling constant 



JV + (0)m 



A {V u + V s ) 

Spp y/V u V s 



(12) 



be- 



For V s /V u <C 1 with fixed Fermi momenta, g c 
comes arbitrarily weak. 

Expanding the gap equation I|1U|) for small A < SpF, 
we find the critical behavior 



(.9 



(13) 



with g c defined in Eq. lfT2")l. 

The quasiparticle excitations of the breached pairing 
phase differ qualitatively from those of BCS. In the BCS 
phase, E p branches are both gapped with their mini- 
mum energies given by |e~ Q ± 2AyJV u V s /(V u + V s )\ ^ 0, 
respectively. By contrast, in the breached pairing phase 
E+ intersects zero energy. Only the E~ branch is fully 
gapped. 

For finite A, the gap equation (|10fl needs to be solved 
sclf-consistently, implementing the restriction in inte- 
gration region, |p| < p^ and |p| > p\, which de- 
pends on A. To this end, we transform the integral 
gap equation into an algebraic one by introducing a 
(BCS) gap paramter, A = A{8p F = 0). The BCS 
gap Aq implicitly defines the interaction strength, l/g = 



3 



JV+(0)ln[A(V u + V S )/A Q }. We shall use A to param- 
eterize the interaction strength below. Performing the 
integrals 



Po+A 



dp 



Po + A 
Po-A 



P "T "0 

dp 



A 2 



pi 



dp 



Pa a/4 2 + A 2 



(14) 



where £+ = ^(V" + V s ) (p — po), we obtain 



A 2 £% 
A = _Pa 

A 2 e s _ 

Pa 



(15) 



for A — > oo. After some elementary algebra, the above 
equation reduces to 



A 2 
A 2 



dp i 



tip%~^ 



tip i 



tip%~^ 



It has a (stable) solution only if 

A 



tipF < 



(16) 



(17) 



which is equivalent to g > g c (see Eq. Ijl2|l ). In solving 
Eq. fTHjl . we distinguish three cases: 

I: Chemical potentials fixed. Then SpF remains un- 
changed before and after pairing, and the gap equation 
permits two solutions: 



Ar 



where Sp 



th 



A y/V"V* (5 PF - S Pth ) 

(threshold value) 



(BCS) 

(unstable BP) 



(18) 

, In this case, 

v v u v s v 

the BP state always has higher energy than the normal 
state, so it is unstable (see Eq. I|24fl h Fig. ^ makes it 
further clear that even if it is a metastable solution at 
tipp > <5pth, the BP state always has a smaller A than 
BCS. It is therefore energetically disfavored. This case 
was considered by Sarma and others |6| a long time ago. 

77: Overall quark density fixed. In this case one has to 
adjust the overall chemical potential, fi oa = \{V u p F + 
V s pp) — \{V U + V s )p , in order to accommodate a fixed 
overall density when pairing occurs. We hold the relative 
chemical potential, <5/i = V u p F — V s p s F fixed. The mis- 
matched Fermi momentum is then altered, SpF — > SpF, 
with '~' denoting the new Fermi momenta p F ' s after pair- 
ing. 

We find Spp is determined by 
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(19) 




25 



20 



15 



10 









(b) 


A Q = 20 MeV / 




A 

BCS / 1 






\Breached 


A =10,MeV 


\Pairing 


A 





50 100 150 

8 p p (MeV) 



50 

5 p p (MeV) 



100 



FIG. 1: Solution of the gap equation as a function of the Fermi 
momentum mismatch Sp F for different coupling constants pa- 
rameterized by Aq . (a) Chemical potentials fixed: The dashed 



line, SpF 



2A 



marks the transition between the BCS 



and breached pairing phases. The range for the unstable 

Ao < 5p F < 



breached phase is ■ 



->A„ 



VV ^ - «v - VV ^- (b) The total den- 
sity fixed: The dashed line, ploted by 5pF = J^, t (1 — a 2 ), 
marks a continuous transition between the BCS and breached 
pairing phases. The range for breached pairing is defined by 

2A yv"V' ) - SpF - vv°v° ■ Both and are P lotted for 
V~ a = 1 and V s = 0.1. 



v^Cylii (tip F - tip F ) ■ For large Fermi 
momenta compared to the mismatch in Fermi surfaces, 



where po = po 



SpF ^PfiPfi we nn< ^ 
SpF = 



tipF 

1 



(20) 



where a = (V u - V S )/(V U + V s ). The gap equation fTfi) 
(with SpF — * tipp) and Eq. I|2(J|) are solved by 



A , 



A 



A VT 7 n^( ^-_^ ) 



(BCS) 
(BP) 



(21) 



The breached pairing phase is realized for 2a 2 — 1 > 
0, i.e., V s /V u < 3 - 2a/2- (See Fig. Ex) The dashed 
line shows a continuous phase transition between BCS 
(tipF _ tipp) and breached pairing (tipF > tipp), where 
Spp = -jLjL= (1 — a 2 ). There is no overlap region where 
both states are solutions. The breached pairing state is 
stable and has lower energy than that of the unpaired 
matter (see $Z%$ below). 

For higher values of V s /V u (but smaller than one), the 
BP solution becomes unstable, being of the same class as 
ifTSjl . Very close to 2a 2 -1 = one finds from Eq. J2U 
that Spp approached <5pth- The value 3 — 2y/2 is model 
dependent, and will be different when a more realistic 
model and a more accurate approximation scheme are 
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used. In the limit V s /V u — > (so a — > 1), the gap 
equation i|16|) exhibits the breached pairing phase for all 
Spp > 0; BCS occurs only at the point Spp = 0. 

In mean field theory there is a single second-order 
phase transition between the BCS and BP phases 
(Fig^>). At fixed mismatch in Fermi momenta Spp and 
as Ao increases, one first encounters BP and then the 
BCS phase. The transition between the BP and BCS 
phases occurs at Spp = 5pp. Since H21|) is symmetric 
with respect to interchange of V u and V s , the expres- 
sion for critical points holds equally for both interior and 
"exterior" gap limits. 

777: Relative quark density fixed. In this case, we also 
hold the overall chemical potential fixed while allowing 
the overall density to vary. In the limit po ^> A ^> A, Spp, 
we find 



4A 2 



(22) 



The gap equation (|16fl (with Spp — > Spp) has no BCS 
but BP solution for any Spp > 0, 



A Q W U V S (S Pth - Spi 



(BP) 



(23) 



The BCS solution only occurs at a single point, Spp = 0. 

Condensation energy. To leading order in A/ Spp, the 
energy difference between the breached pairing and nor- 
mal (A = 0) phases is 



Erp — Em = < — 



W + (0)A 4 
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2Spj, 


yuys , 


W+fOJA 4 
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(I) 
(II) 
(III) 



(24) 



The BP state is energetically favorable in case II if 
V s /V u < 3 - 2V2 (with overall density and Sp, fixed), 
and it is always favorable in case III (with the relative 
density and overall \i fixed). 

In the breached pairing scenario pairing occurs at zero 
total mometum and peaks near the light species Fermi 
surface. It resembles standard BCS in that the pairing 
occurs at all angles, there is no necessity for translation 
or rotation symmetry breaking, and there is a branch of 
gapped excitations. There are, nevertheless, also gapless 
modes at two new "effective" Fermi surfaces, between 
which pairing is suppressed. For small A the energy gain 
of the breached pairing phase with respect to unpaired 
normal matter is ~ A 4 , parametrically less than ~ A 2 in 
the BCS case. 

Discussion. In the context of QCD, enforcing charge 
neutrality will make case III relevant for the u and d 
quarks. For pairing between strange and light quarks, 
case II is relevant. Thus, breached pairing allows all 
flavors to form condensates (ud) bp, (us) bp, (ds) bp of 
breached pairing type. We suspect that there is a density 



range where it is energetically favorable compared to un- 
paired matter or the 2SC phase with only (ud)scs 7^ 
condensation. Alford and Rajagopal Q showed the ab- 
sence of 2SC phase in compact stars, and suggested that 
there should be some non-BCS pattern of pairing includ- 
ing possible unpaired matter in the range of densities 
between hadronic and CFL matter. The breached pair- 
ing version of CFL appears to be a serious candidate for 
these intermediate densities. Since there is a simple sec- 
ond order phase transition from BP to BCS superfluidity, 
as a function of increasing density, it is natural to expect 
a single phase transition from breached pairing CFL to 
standard CFL quark matter. Detailed calculations, in- 
cluding shifts of Fermi surfaces due to electric and color 
neutrality, will be required to determine the range, if any, 
over which breached pairing CFL is favored (in prepara- 
tion). At strong coupling one could even speculate about 
breached pairing CFL states at zero strangeness, where 
the strange chemical potential falls between ±m s . 

The possibility that one might rearrange the free- 
particle Fermi surfaces to take advantage of opportunities 
for pairing was considered long ago in condensed matter 
literature by Sarma and others [6j. Superfluidity with 
quasiparticle dispersion relation crossing zero, thus lead- 
ing to gapless modes, was encountered in a model re- 
sembling ours in Ref. In neither case, however, was 
the possibility of a stable breached pairing state recog- 
nized. (See also a recent study in cold ato ms H .) As this 
work was being completed reference Ref. appeared, 
in which a related problem, imposing a charge neutrality 
condition directly, was analyzed. 

We have phrased our discussion in terms of QCD and 
quarks, but one can easily imagine closely related univer- 
sality classes for condensed matter or cold atom systems. 
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